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Abstract 

It has been argued by Grigoriev and Rubakov that one can simulate real time 
processes involving baryon number non-conservation at high temperature using real 
time evolution of classical equations, and summing over initial conditions with a 
classical thermal weight. It is known that such a naive algorithm is plagued by 
ultraviolet divergences. In quantum theory the divergences are regularized, but the 
corresponding graphs involve the contributions from the hard momentum region and 
also the new scale ~ gT comes into play. We propose a modified algorithm which 
involves solving the classical equations of motion for the effective hard thermal loop 
Hamiltonian with an ultraviolet cutoff fi 3> gT and integrating over initial conditions 
with a proper thermal weight. Such an algorithm should provide a determination of 
the infrared behavior of real time correlation function < Q{t)Q(0) >t determining 
the baryon violation rate. Hopefully, the results obtained in this modified algorithm 
would be cutoff-independent. 
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1 Introduction 



Sphaleron processes in electroweak theory mediate baryon number violating pro- 
cesses, and have been computed to be significant at high temperatures [Q — [|]. 
This computation involves estimating the probability of a classical thermal fluctua- 
tion evolving over the top of an energy barrier which separates topologically distinct 
phases of electroweak theory. At very high temperatures, where the masses of the 
W and Z boson can be ignored, the field configurations which dominate the contri- 
bution have a typical size of order l/oty/T |J . This is the length scale on which 
non-perturbative physics is important, and there is no known way of doing a weak 
coupling estimate of the effect. 

In a provocative paper by Grigoriev and Rubakov , a numerical algorithm was 
suggested for the estimation of the baryon number violation rate. They argue that 
at high temperature, the dominant process for sphaleron processes involves transi- 
tions of fields over the energy barrier. The dominant transitions should therefore 
be classically allowed thermal processes. Since the field configurations useful for 
sphaleron processes are long wavelength and involve many quanta, it is reasonable 
to assume that the evolution of the field configuration over the top of the barrier is 
in fact classical. 

If the evolution of the fields is classical, then they will be described by classical 
equations of motion with given initial conditions. These initial conditions are de- 
termined by a classical thermal ensemble where the initial distribution of fields and 
their conjugate momenta are computed in a Monte-Carlo simulation. 

This method was used to compute the rate of sphaleron transitions in elec- 
troweak theory at temperature higher than the electroweak phase transition tem- 



perature ||. It was also tested in two dimensional models || — [fL0|| . Although 
the two dimensional results looked convincing, there have remained a number of 
plaguing problems associated with obtaining correct dependences on effective tem- 
perature dependent masses, and these masses depend on the ultraviolet region of the 
theory where, as we will see, there are problems associated with treating the theory 
classically. The 3+1 dimensional simulations of electroweak theory are plagued by 
ultraviolet singularities associated with making the spatial lattice spacing smaller, 
and one may worry whether any quantitative numerical conclusions are warranted. 

It is easy to understand why there are problems associated with taking the zero 
lattice spacing limit in such simulations. In the second section, we shall present a 
detailed mathematical argument in the context of theory . In the remainder of 
the Introduction, we shall present heuristic arguments which display the nature of 
the problem. 

It has been well known since Rayleigh, Jeans, and Planck that, if one tries to 
describe a bath of photons classically, the free energy is ultraviolet divergent and 
depends on the cutoff. To see this, consider the energy density for a gas of scalar 
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massless bosons. This density at high temperatures behaves as 



(2tt) 3 e^IPl - 1 30 

On the other hand, if we try to compute the thermodynamic characteristics of the 
system classically, we have the functional integral representation 

Z d = J [d<j)}[dn] exp{-(3H} (1.2) 

where the Hamiltonian is 

H = J Sx{hi{xf+ l -{V<P{x)f] (1.3) 
The results are ultraviolet divergent. For the energy density we get 

E d (T) = -^/i 3 (1.4) 

where \i is a momentum ultraviolet cutoff (The simplest way to get this result is to 
substitute the Bose distribution function in Eq. (|1 . 1| ) by its classical limit T/|p| ). 
This is the well known classical Rayleigh- Jeans divergence and was one of the original 
motivations for quantum mechanics. 

Ultraviolet divergences also crop up in other quantities. Consider e.g. the 
correlator < <j>(x, t) 0(0) >t in 4 theory. 

This bears on what we discuss here since the quantity of interest defining baryon 
nonconservation rate in the standard model is a similar ( but, of course, more 
complicated) correlator 

C(t) =< [Q(t) - Q(0)]Q(0) > T (1.5) 

where 

Q{t) = gf^ey* / d 3 x (F%Al - 9 -e abc A1A)Alj (1.6) 

is the Chern-Simons winding number of the SU(2) gauge field configuration. The 
expectation is that the correlator ( [[ . 5|) behaves as C(t) ~ FV\t\ at large \t\ 3> 
(awT)^ 1 where T is the baryon number nonconservation rate per unit volume and 
is estimated as 

r A ^°(T) = n{a w Tf (1.7) 
We take the interaction term in the action as 

S = J dt J d 3 x ^0 4 (x) (1.8) 
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This yields in lowest order, the diagram shown in Figure [1|. 

In this paper, all perturbative calculations are performed consistently in the 
real time formalism which is the most suitable tool to study kinetic properties in 



thermal medium. In its complete form the formalism has been developed in [11 



I2fl . It involves some intricacies, one should carefully distinguish the retarded, the 
advanced and the so called P-component of the Green's function []. But we do not 
really need to go into details here because, for our purposes (we will be interested 
only with 1-loop graphs and only with the real part thereof), it almost suffices to 
use the simplified version of the formalism due to Dolan and Jackiw [|17] (what 
almost means is explained in the Appendix A). In Dolan-Jackiw formalism the 
thermal propagators present the sum of two terms: the standard zero-temperature 
propagator and the temperature-dependent term which reflects the presence of real 
particles in the heat bath (alias, the temperature insertion). In particular, the 
thermal propagator of massless scalar is 

D T (p) = ^— + 2^)^—- (1.9) 

The temperature-dependent contribution in the polarization operator in Fig. [I] 

is 

X ' d3 » ^ (1.10) 



2 J (27r) 3 |p|(e^lpl - 1) 24 

which is the familiar temperature renormalization of the scalar mass. (Here and in 
the following, we ignore the infinite ultraviolet mass renormalization in the zero- 
temperature theory which we assume to be cancelled by a proper counterterm.) Let 
us try to calculate the same graph in classical theory. It is ultraviolet divergent Q. 



Note that the method of 0-|TIJ where the correlator ( |1.5| ) was estimated by solving 
the classical equations of motion with the given initial conditions which then were 
averaged over by integrating over classical phase space with the weight exp{— /3H d } 

The review can be found in H, @. For applications of the formalism to hot gauge theories 

see e.g. @ |H§ . 
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Up to now, we assumed that the tree- level mass of the scalar field is zero, but the result (1.11) 



holds also for the theory with nonzero mass Mq if /u, !2> Mo- For the one-loop estimate (1.11) to be 
stable with respect to higher-order corrections, the characteristic momenta in the integral |p| ~ fi 
should also be much larger than the full quantum temperature-induced mass 
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like in ( |1.2| ) is equivalent in perturbation theory to calculating the classical graphs. 
In Sect. 2, we shall illustrate this explicitly for A0 4 theory. 

The way out of this difficulty suggested in |J was to choose /i ~ T. The 
particular proportionality coefficient was fixed by the requirement that the classical 
and quantum expressions for the energy density (the analogs of (|1.1| ) and ( |1.4|) ) 
coincide. Note, however, that choosing /i in such a way, we would still have a 
mismatch between the classical ( 1 . 1 1| ) and quantum ( 1.10|) expressions for the mass 



renormalization. This mismatch is a manifestation of a simple fact that the classical 
description is justified only for the low momentum modes. But, as we have seen, the 
low-energy sector of the theory does not entirely decouple from the high momentum 
modes. They get entangled together. Thus, the procedure used in || may be not 
reliable []. 

Our main observation, however, is that low momentum modes and high mo- 
mentum modes can in principle be disentangled. The correct procedure should be 
the following. Let us introduce an intermediate scale /i such that gT « /i < T 
and treat the modes with momenta |p| > /x and |p| < fi in a different way. The 
high momentum modes can be explicitly integrated over leading to the effective 
hard thermal loop Hamiltonian for the low momentum modes As the low 

energy bound for momenta circulating in the loop is not zero but is equal to fi, 
this effective Hamiltonian would involve also some counterterms depending on the 
parameter fi. With the effective Hamiltonian in hand, one can return back to the 
original GR procedure, solve the effective equations of motion, and average over ini- 
tial conditions with the weight exp{— (3 H e " (phase space)}. The ultraviolet 3-dim 
momentum cutoff should be chosen to be of order of /i. The final result will be 
/i-independent []. 

We emphasize that the success of this modified procedure depends on two cir- 
cumstances: 

1. For high momentum modes, the classical treatment is wrong, but these modes 
can be taken into account perturbatively. 



3 This question is, however, not clear at present. We return to the discussion of this issue in the 
last section. We want also to emphasize that the violent power ultraviolet divergences in m cl (T) 
and other parameters of the effective hight-T Hamiltonian are spesific for 4-dimcnsional theories. 
In two-dimensional case considered in the first Grigoriev and Rubakov paper, the situation is much 
more benign. 

4 We were able to construct such a /^-independent algorithm explicitly only for the scalar theory. 
In Ref. [^lj such an algorithm has been discussed for a scalar theory in (2+1) dimensions. The 
implementation of this idea in the gauge theory meets serious technical difficulties. We shall discuss 
it at length later. 
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2. For low momentum modes, perturbation theory breaks down and the classical 
real time equations of motion are essentially non-perturbative. They cannot 
be evaluated by weak coupling methods. But the classical approach is justified 
here, and one can do a numerical calculation in the GR spirit. 

Note that the method of separating hard and soft modes where the former can 
be treated perturbatively while the latter cannot, but the influence of hard modes 
on the soft mode physics is taken into account, is a standard one in many physical 
problems. E.g., it is the basis of the operator product expansion technique and sum 
rules in QCD [|| . 



In this paper, we first restrict ourselves with an illustrative analysis for A0 4 
theory. Later we outline the procedure to be used for gauge theories. Note that 
the situation is much more complicated in gauge theories when one computes the 
real time evolution of classical fields. For computing the real time evolution of the 
gauge fields, one obtains Feynman diagrams of the type shown in Fig. ||]. When 
these diagrams are computed, one obtains a correction of the form 

aT 2 F^(k°/\k\) (1.12) 

where /c M is the external momentum. This function has a nontrivial dependence 
on A;°/|k|. For |k| ~ gT, this contribution is as large as the free particle kinetic 
energy terms, g^ u k 2 — k^k v . One can show, however, that so long as k° > gT 
and |k| > gT the contribution of higher-order loops in the two-point function and 
other irreducible vertices is still suppressed and perturbation theory works. The 
shorthand for all these vertices is called the Braaten-Pisarski effective Lagrangian 



IS 



The main problem which distinguishes hot gauge theories compared to the sim- 
ple A0 4 case is that this effective Lagrangian is highly nonlocal. One cannot straight- 
forwardly write the effective theory in a Hamiltonian form and write differential 
equations of motion for real time evolution of initial thermal gauge field configu- 
rations. However, in the recent works [|T^, |2D, ^] this theory has been written in 



the local Hamiltonian form. The price which one had to pay is the introduction 
of additional variables the integration over which restores the original nonlocalities. 
Still, the local Hamiltonian description of hot gauge theories exists. Unfortunately, 
such a local effective Hamiltonian exists only if performing the integration in the 
hard thermal loops over all loop momenta (without infrared cutoff fi). As was al- 
ready mentioned, we were not able to derive /^-dependent piece in the effective HTL 
Hamiltonian. Our hope is, however, that the contribution of these unknown terms 
in the baryon nonconservation rate in interest is suppressed. 

The organization of this paper is as follows: In the second section, we reconsider 
the formalism of Grigoriev and Rubakov, and show explicitly how the ultraviolet 
cutoff dependence arises. In the Section 3 we use this formalism to compute the 
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lowest order contribution to the propagator in the A0 4 theory. In the fourth section 
we present a path integral formalism for computing real time correlation functions 
by a method which in lowest order approximation gives a result similar to that of 
Grigoriev and Rubakov. The effects due to hard modes are discussed in Sect. 4. We 
argue that, for most practical purposes, one can evaluate the contribution of high 
momentum modes in a Gaussian approximation to the functional integral represen- 
tation for the real time correlation function. We explicitly demonstrate this for \<j) 4 
theory. We show that this results in a formalism similar to that of Grigoriev and 
Rubakov except that one should use the effective hard thermal loop Lagrangian plus 
some counterterms needed to render the three dimensional classical equations finite. 
In Sect. 5, we discuss hot gauge theories and, in particular, the kinetic equation ap- 
proach and the Nair's effective Hamiltonian which may be the basis for the correct 
analysis of the problem of interest. We also discuss the kinematic limits in which 
our computation is valid. In the last section, we summarize our results. We describe 
what will be needed to be done before realistic computations for 3+1 dimensional 
gauge theories may be performed. 

There are also two technical Appendices devoted to two unsuccessful attempts 
to isolate /z-dependent terms in the effective Hamiltonian for hot Yang-Mills theory. 
In Appendix A we isolate the leading cutoff-dependent terms in the soft multigluon 
thermal vertices with momentum infrared cutoff. Unfortunately, such a cutoff breaks 
gauge invariance. As a result, the fi - dependent terms in the effective lagrangian 
are not gauge invariant, and the classical problem cannot be consistently posed. In 
Appendix B we display the problems one meets when trying to do the same with a 
lattice cutoff (here the problems are of the opposite kind. The lattice regularization 
is gauge-invariant, but to isolate the /z-dependent pieces is very difficult technically) 
and perform some illustrative calcultations for hot scalar QED. 

2 Ultraviolet Problems 

In this section, we elucidate why the procedure of Ambjorn et al. doesn't properly 
take into account the high momentum modes with k ~ T. This is demonstrated for 
the case of a scalar theory with the action 




(2.1) 



The object we want to compute is the real time correlation function 

C(*,k,p) =< 0(t,k)0(O,p) > T 



(2.2) 



where 




(2.3) 
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In general, one can consider correlation functions of composite operators, but for 
the points we wish to make, it is sufficient to consider correlation functions of the 
elementary boson field. 

The procedure of Grigoriev and Rubakov is equivalent to the following: Assume 
an initial condition which has been determined by random selection from a classical 
thermal distribution, that is, with weight 

exp{-(3H(U h \r)} (2.4) 

where II is the momentum conjugate to 0. If we know Il m and m , we have the 
initial conditions needed for solving the classical equations of motion. 

The next part of the recipe is to solve the classical equations of motion. This 
then determines U(t, k) and 0(t, k). To compute the correlation function, one then 
has to average over the initial conditions. 

The equation of motion 

(^ + k 2 + M 2 ) 0(t,k) 

= "fl / (2^ (2^ <W> k ^ k ^ k " kl " k2 ) ( 2 - 5 ) 
is solved perturbatively: 

0(t,k) = 1 (t,k) + 2 (t,k) + --- (2.6) 
0i (t, k) is a superposition of plane waves: 

0x(t, k) = — (a(k)e-^ fc '- kx) + a*(-k)e^'- kx) ) (2.7) 

where a, a* are determined by the initial conditions and 

oo k = y/k 2 + Ml (2.8) 
The mass M c \ is related to the tree level mass by 

Ml = Ml + 5Ml (2.9) 
The mass counterterm 5 Ml will be determined below. 02 (t, k) is obtained by solving 

g + k 2 + Ml)0 2 (t,k) = -y^^M^MM) 

x (t, k - k x - k 2 ) + 5MlMt, k) (2.10) 
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with the initial conditions 

2 (O,k) = O, ^> 2 (0,k)=0 (2.11) 

To evaluate the Green's function (|2.2|), we have to average m (p) [<fii(t, k) + </> 2 (£, k)] 
over the initial condition with the weight ( |2.4fj : 

c d (t, k, P ) = z c f / [rfn in ] [rf0 in ] exp{-/3i/ [n in , <T]} 



where 



^4! v 7 - / 2"" cl 
(p)U(t,k) + 2 (t,k)) (2.12) 



F o [n,0] = / rf 3 a;Un 2 (x) + ^(x)(-V 2 + M c 2 1 )0(x)| (2.13) 



One obtains 



C d (t,k,p) = (27r) 3 5(p + k)^|cos(w p t) 



~2^{ XT I 1^~ 6M *) K tsin (M + cosKt))|(2.14) 
The momentum integral is ultraviolet singular and is regularized by the cutoff yu: 

« iL^M, (2 , 5) 

The contribution proportional to t on the r.h.s. of eq. ( 2.14Q is due to a resonance 
term in the r.h.s. of eq. ( |2.10|) and has to vanish. This determines the mass countert- 
erm (this self-consistent procedure is just equivalent to solving the Dyson equation 
with the polarization operator calculated by the tadpole graph in Fig. [I] in the 
classical limit): 

^ci=^ + 0(-|l) (2 . 16) 

Taking the Fourier transform of C c \ with respect to t we obtain the classical Greens 
function in momentum space as 

G cl (p) = j^2rr6(p 2 - M 2 ) (2.17) 
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This has to be compared with the quantum mechanical expression. The summation 
of the diagrams in Fig. ^| yields the real time propagator 



with 

M 2 = M 2 + 

24 



M 2 = Mq + — T 2 (2.19) 



The low energy limit eq. (|2.18 



is 

G(p) ~ ^2tt5(p 2 - M 2 ) (2.20) 

(It is also the classical limit. When substituting p — > hu, p — > ^.k, we see that the 
second term in Eq. ( |2.18| ) involves an extra power of h in the denominator compared 
to the first term when — > 0). 

Ambjorn et al. have choosen the cutoff such that the energy densities of the 
quantum mechanical (eq. fll.lQ ) and the classical system (eq. ( |1.4|) ) are equal. With 
this choice, however, there is a mismatch of the self energies corresponding to the 
mass counterterms. 



3 A Formal Solution 

As was emphasized in the introduction, the quantities of interest which we wish to 
compute are real time correlation functions 

C(t,x) = (0(t,x)0(O)) T 

Tr{e-P H e itH O(0, x)e~ UH O(0)} 



Tr e~P H 
Tr{e im e- pH O(0, x)e"^C»(0)} 



(3.1) 



Tr e~P H 

where O[0>(t, x)] is a local operator. This correlation function can be expressed as 
a path integral as 

C{t,x)= f [dll\[d<5>} O{t,x)O(0)e iS (3.2) 

• / *Uo = _ i/3 =*| a .O =0 

In this equation, the time variable t is a real Minkowski time. The fields are periodic 
over the time x° = to time x° = —i/3. The action has been written in terms of 
fields and their conjugate momenta in the form 

S = J^dx J d 3 x (n(x)d $(x) - H{<$>, V$)) (3.3) 
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where Ji is the Hamiltonian density. The contour C of integration in time is shown 
in Figure |j. It begins with a real time evolution from x° = to the time x° = t. 
It then takes a Euclidean path from x° = t to x° = t — i(3. Finally x° evolves 
backwards in time from x° = t — if3 to x° = —i(3. The periodic boundary conditions 
for the trace are at the end points of the time contour. 

This path integral representation yields the expression of Grigoriev and Rubakov 
in the stationary phase approximation 5lm(iS) = 0. This requires IT and $ to 
be constant along the Euclidean interval of the contour, i.e. n(x°,x) = IT m (x), 
$(x°, x) = $ m (x) and to be a solution of the classical equations of motion along the 
real time pieces. Then the phase vanishes in region Cu, while the phases of regions 
Cj and Cju cancel each other. The exponential exp(iS) becomes exp{— /3H(H m , $ m )} 
and the only remaining integrations are those over IT 11 and <3> m , which are the ini- 
tial values for the real time evolution in regions Cj and Cm. Then the correlation 
function is 

C d {t,x) = J [rfn in ][rf$ in ]e-^( nin '* in )O d (t,x)aK0) (3.4) 

This treatment ignores quantum corrections completely. It can be justified if the 
only relevant modes are those with momenta much smaller than the temperature 
since then their occupation number is large and they can be treated classically. 
Since the calculation of Grigoriev and Rubakov is performed on a lattice, it involves 
a cutoff /i ~ 1/a, where a is the lattice spacing. From the above argument fi should 
be much smaller than the temperature. Grigoriev and Rubakov have choosen fi ~ T. 
If their result is //-independent, it is probably reliable. However, if it depends on /i, 
it is sensitive to modes with momenta of order T, for which the classical treatment 
is not correct. 

For the latter case we suggest the following procedure: One should choose the 
cutoff /i such that 

gT<^fi^:T (3.5) 

(This condition is written for gauge theories. For the A0 4 theory the characteris- 
tic temperature-induced mass gT should be substituted by A^ 2 T). The fields are 
decomposed into short wavelength modes and long wavelength modes: 

n = n 5 + n L , $ = $ 5 + $l (3.6) 

where ITs and $5 contain the Fourier componets with |k| > /i, and 11^, $l consist 
of those with |k| < /i. Then the $5 interact weakly among themselves and they can 
be treated as free fields in the background of i.e. for these modes one can use 
the Gaussian approximation in the path integral. On the other hand, the condition 
ji < T ensures that the modes with |k| < /x can be treated classically. 
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Naively we could ignore the Gaussian fluctuations around a classical solution 
because it is higher order in coupling and, by assumption, the coupling is weak. 
In this case, however, there are corrections of order g 2 T 2 to the classical equations 
of motion and these generate a huge correction for the small momenta we wish to 
consider. These terms arise from momenta which are much larger than the typical 
frequencies we wish to consider for our classical solution, and therefore the Gaussian 
fluctuations are dominated by the ultraviolet. The integration over these modes is 
equivalent to integrating out the high momentum modes and generating an effective 
theory of the low momentum modes. 

Let us examine in more detail the reason why it is a good approximation to 
only do the Gaussian corrections to the path integral in the presence of the long 
wavelength modes. Consider modes with a typical momentum scale //. In this range 
of momentum, the occupation number of the modes, T/fj,, is large compared to 1, 
but small compared to 1/g. When these modes interact among themselves, the 
typical strength of interaction is of order gT/^ <C 1. Thus, nonlinear interactions 
among hard modes can be neglected, and the latter can be treated in the quadratic 
(Gaussian) approximation. In a diagrammatic language, it suffices to keep track 
only of the graphs involving 1-loop subgraphs with high internal momenta (the so 
called hard thermal loops). 

We therefore seek to integrate out the short wavelength modes to obtain an 
effectively classical theory for the long wavelength modes. This gives 

C(t,x) ~ Acm L ][d$ L ]e i(s[ni '* il+55[ni '* i; ^C>[$ i (t,x)]C)[$ I/ (0)] (3.7) 

where 

e i5S{n L ,4> L ;ti] 

= J [du s ][d^ s ] exp ji Jdx° (J d 3 x (n s d $s - \n 2 s - 5H{ ^ l) ^(x)^ 

~ 5/^/^ w "^g y) *^.W,y)))} (3-8) 



The effective theory described by Eqs. (|3.7| , ^7§|) will be the one for which the GR 
treament will be applied. 



4 Integrating Out the Hard Modes 

We have seen in the last section, that the effect of the hard modes is to generate 
an effective action for the long wavelength modes. Furthermore, it is sufficient to 
integrate them out in the Gaussian approximation. On the sight, this still appears as 
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a tremendously difficult task. The exponent on the r.h.s of Eq. ( |3.8| ) contains terms 
linear in They can be nonvanishing when contains modes with momenta 
smaller but close to /i and the fields would have to be shifted, $5 — > $5 + 8$s- 
Furthermore one has to sum all one loop diagrams with external momenta |pj| < //. 
The loop integration over k involves theta functions of the type 0(|k + Pj| — /i), 
where the Pj are linear combinations of the external momenta, for each internal 
line. However, these complications arise only for relatively hard modes in $^ with 
I Pi I ~ \i. In the previous section we have argued, that these modes are only weakly 
interacting and the corresponding contribution to SS should be small. f\ 

Therefore we can restrict ourselves to external momenta \pi\ ~ gT <C //. In 
this approximation SS\)Jl, $l] A*] is given by the sum of the one loop diagrams with 
a factor <&l{x) for each external line, where the loop momentum is restricted by 
|k| > /i. Then it is possible, even for a nonabelian gauge theory, to extract the 
leading terms of all n-point functions. This is similar to the analysis by Braaten 
and Pisarski |Tj|, who extracted the leading contribution proportional to g 2 T 2 . The 
only difference is that in our case there is a lower limit /1 for the |k| integration, 
where k is the loop momentum. 

For illustration, let us consider the simple case of the A0 4 -theory with the action 



1 



c J w rw 2 



S= / dx° / d 3 x { 7i{x)d (f){x) - - 7T 2 (x) + (V(j)(x)) 2 + m 2 <J) 2 (x) --7<j> 4 (x) (4.1) 



A 



4! 



Here the leading contribution in the effective action is given only by the diagram in 
Fig. |l[ First the integration over tt(x) is performed. Then the free propagator for 
the short wavelength field is 

iDs(z 1 ,x i ) = J ^e(|k|- M )po(A;)e-^-^(e c (x?-x°)+n(fc )) (4.2) 

where n(k ) denotes the Bose distribution function, Oc(^? — ^2) i s ^ ne s ^ e P function 
along the contour, and the spectral density po(k) is given by 

Po(k) = 2n^- (6(k° - u k ) - 5(k° + u k )) (4.3) 

with u k = (k 2 + m 2 ) 1 / 2 . We obtain 

A 



6S=-i , 
2 Jc 



dx° J d 3 x(p 2 {x)Ds{x,x) (4.4) 



In other words, the exact propagators and vertices with external momenta which are much 
larger than the characteristic temperature-induced mass coincide with tree propagators and vertices 
plus small corrections. 
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The value of Ds{x, x) is independent of x°. In particular, it is the same on all three 
pieces of the contour. Therefore Eq. (4.4) can be written as 

5S = f dx° J d 3 x (--8m 2 (j) 2 (x)^ (4.5) 

This expression contains an ultraviolet singularity, which is removed by a mass 
renormalization at zero temperature. Then, in the limit m 2 <C « 2 , one obtains 

Sm 2 = A (t 2 - A/iT) (4.6) 



24 V vr 

Taking into account Eq. (|4.5| ) when calculating the correlation function has two 
effects. First, it modifies the equation of motion for the real time evolution in 
regions Cj and Cm of the time contour. Secondly, it changes the Hamiltonian on 
the Euclidean piece Cjj, which gives the thermal weight when taking the average 
over the initial conditions. 

For A0 4 theory, there are two equivalent ways to impose an ultraviolet cutoff in 
the classical problem. First, one can just expand the field in a Fourier series 

0(x,t) = $> n (t)e 2 ™~ l\ (4.7) 

n 

write the classical equations of motion in terms of Fourier harmonics n (t) and 
impose the momentum cutoff |k| = 27r|n|/L < «. Technically, it is more convenient, 
however, to define the theory on a Hamiltonian lattice. As we show in Appendix B, 
the ultraviolet divergences in the one- loop lattice graph can be easily singled out, 
and the lattice analog of Eq. ( |1.11| ) is given by 

2 /rrn ATa 2 f 7T / a dpx f*/ a dp 2 f' K l a dp 3 1 



M 2 cLlat {T) 



<r/a 2lT J-ir/a 2lT J-n/a 27T 3 — Z)f=l COs(apj) 

XT f n r n f 77 I XT r°° 

dx dy dz- = irl ^« e" 3 % 3 (a) (4.8) 



4(2ir) 3 a J-tt J-n J-n 3 — cos x — cosy — cos z 4a 

where a is the 3-dimensional lattice spacing and Io{ct) is the exponentially rising 
variety of Bessel function. We see that the proper choice of a is 



2 

da e- 3 % 3 (a) (4.9) 



TT 

H Jo 



Then the formulae ( |4.8|) and ( 1.1 1|) match each other and the final result for the 
correlators fl3.1|) is cutoff-independent. 
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5 Gauge theories 



Let us discuss now whether it is possible, and if yes then how, to generalize the 
method worked out in the previous sections for A0 4 -theory to an SU(iV) gauge 
theory. [] As has already been mentioned in the Introduction, this is not so trivial 
since integrating out the modes with |k| ~ T does not simply result in a mass 



renormalization as in Eq. fll.10 ). One rather obtains nonlocal terms proportional 
to g 2 T 2 , which for the modes with |k| < gT are at least as large as the terms in 
the classical action (for the modes with |k| ~ g 2 T they are even much larger!). 
Furthermore, there is a leading g 2 T 2 contribution to all n-point functions. These 



terms have been identified in Ref. |T3] and their generating functional is known as 



the hard thermal loop effective action. Since the resulting equations of motion are 
nonlocal, they are not appropriate for computing correlation functions a la Grigoriev 
and Rubakov. 

Fortunately, this hard thermal loop effective field theory can be presented in a 



local form [|T9[ |2(J . In Ref. [19[] it was formulated as a set of kinetic equations for 
the long wavelength gauge fields and the induced current j of hard particles moving 
in the slowly varying background: 

D,F^ = f (5.1) 

f(x) = 2gN J -0-v»5n(p,x) (5.2) 

Here denotes the covariant derivative in the adjoint representation, 5n = 
5n a t a is the deviation of the thermal distribution of the hard particles from the 
equilibrium Bose distribution function no(|p|). Furthermore, is the velocity of 
massless particles with momentum p, = (1, v), v = p/|p|. The kinetic equation 
for Sn is f\ 

v»D,Sn(x, p) = -gv - (5.3) 



5n depends nontrivially only on v and one can integrate eq. (|5.3|) over |p|. On the 



3 We shall ignore for time being the Higgs fields and fermions. Possible effects due to their 



presence will be discussed at the end of the paper. Note that the studying of the correlator (|l.5| ) 



makes sense also in hot QCD where it determines the rate of processes with chirality violation. 
7 The equation system - ) presents the nonabelian analog of the well-known Vlasov equa- 



tions derived originally for the usual non-relativistic plasma |13| . Silin generalized them for ultra- 
relativistic abelian plasma p4| . In Ref.|^] the kinetic equation approach has been used to study 
the properties of supersymmetric field theories at finite temperature. 
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r.h.s. of equation ( |5.3| ) one has the integral 



, 2 dno(|p|) . 1 2 _2 



Then the induced current can be written in terms of 

3 



V ' ' £7T 2 T 2 

as 



J 8n(x, p)|p| 2 <i|p| (5.5) 



f(x) = 2m 2 (T) J -^-v u w(x, v) (5.6) 

where m 2 (T) = |u;^(T) , and ui^iT) = g 2 T 2 N/9 is the plasmon frequency. The 
function w satisfies the equation 

v^D^wix, v) = v • E(x) (5.7) 

By solving eq. ( |5.7|) the induced current can be written as a functional of the 
gauge fields. Then the eq. ( |5.1| ) takes its nonlocal form obtained from the hard 
thermal loop effective action. The energy of the system of gauge fields and induced 
currents is given by 

H = J d 3 xTr |e ■ E + B • B + 2m 2 (T) J ^-w(x, v)w(x, v) | (5.8) 

An alternative local formulation has been given in Ref. pO] by introducing aux- 
iliary fields G(x, v), which are SU(iV) matrices satisfying the constraint G(x, — v) = 
Gfi(x, v). In the gauge = with D± = d ± v ■ D, the Hamilton equations of 
motion are: 

E? = d A« 
doEt + evkiDjBtf = JdWiiJl-Jl) 

m 2 (T) 

(D-J + ) a = -— { -^-E« Vi (5.9) 

07T 



where 



J« = H^lT T { lt a (D + G)G- 1 } 

ja = _Hmi T ^ lt a G -l^ D G ^ ( 51Q ) 

2n 
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The Gauss law constraint 



(DE 



\ a + J dv{Jl + r) = (5.11) 



should be imposed. The two local formulations of the hard thermal loop effective 
theory given by eqs. ( |5.1|JT6| , |5T^ ) and eqs. (^9] - |5.11| ) have been shown to be 



equivalent in Ref. [23 



An important remark is in order here. The one-loop graphs provide the dom- 
inant contribution to the effective Lagrangian when the external momenta are of 
order gT. When the external momenta become as small as g 2 T (and we presumably 
need that to solve the problem), the contribution of certain higher-order graphs to 
the effective vertices is of the same order as the contribution of hard thermal loops 
of Figs. |J pi. 

Really, consider the two-loop graph for the gluon polarization operator in the 
kinematic region where the external momentum q is of order g 2 T (see Fig. [5|). The 
extra loop brings about the small factor g 2 , but one can show that this loop involves 
the "resounding denominators": the propagators with momenta p and p + q bring 
about the large factor 

T A 5(p 2 )-^- 2 ~-~- 2 (5.12) 
(p + q) 2 pq g 2 

which compensates the smallness ~ g 2 so that the contribution of the graph in Fig. |5] 
is unsuppressed. (An accurate analysis shows that such a large contribution appears 
only if taking into account hard thermal loops in the internal gluon propagator.) 

In fact, all higher- loop ladder graphs contribute on the same footing. To find out 
the gluon polarization operator, the 3-point, 4-point etc. vertices at small momenta 



~ g T, one has to resum all these ladders which is a formidable task (in Ref. [|15 
such a resummation has been done for the photon polarization operator in abelian 
theory) . 

Fortunately, however, we do not need to do that. The matter is that the integral 
for the graph in Fig. B is saturated by the region where the momentum of the rung 



of the ladder is soft: |k| ~ gT JT^|. In our procedure, the integration should be done 



with the cutoff |k| > [i ^> gT for all loop momenta. In that case, the contribution 
of the graph in Fig. [| and all its higher-loop counterparts is suppressed again and 
can be safely ignored. Thus, we are left with the standard one-loop graphs in Figs. 
0, [7] with loop momenta being restricted to be greater than the separation scale \x. 
Our suggestion is to solve numerically the equation system (|5.9| - |5.1lD rather 



than just equations of motion of the free Yang-Mills Hamiltonian. This would pro- 
vide a proper account of the potentially dangerous high-momentum contributions. 
However, if one just uses the effective Hamiltonian ( |5.8| ) obtained after integration 
over all momenta in the hard thermal loops, the results would, generally, depend 
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on the ultraviolet cutoff which has to be imposed. If one wants to construct a 
cutoff-independent algorithm, one should add to the Hamiltonian ( |5.8|) \i- depen- 
dent counterterms so that the final result would not depend on /i. However, to find 
out these counterterms explicitly is much more difficult task here than in the toy 
A0 4 model. 

For example, a simple-minded momentum ultraviolet cutoff is not suitable as 
it breaks gauge invariance. If trying to include the cutoff-dependent terms in the 
effective Hamiltonian (qualitatively, they modify the thermal mass m 2 (T) ~ g 2 T 2 
by the terms ~ g 2 fiT), such a Hamiltonian would not commute with the Gauss law 
constraint and the classical problem just cannot be posed. See Appendix A for the 
detailed discussion. 

Also, we do not know how to do it with a gauge-invariant lattice ultraviolet 
cutoff. First, it is not so trivial to construct a lattice version of the Nair system 
( |5.9[ - |5.1lD . But, even assuming it might be done, the major problem is that the 
lattice regularization breaks rotational invariance. As a result, the counterterms 
are also not rotationally invariant (see Appendix B for an illustrative calculation in 
scalar QED), their structure is complicated, and we are not able to write down their 
general form. 

Thus, the only practical recipe we can suggest is to take the effective equations 
of motion for the full HTL Hamiltonian ( |5.8|) and solve them numerically in the GR 
spirit with an ultraviolet cutoff fi ^> gT. After averaging over initial conditions, 
we would get the correlator (|1.5|) .Though the result would generally depend on 
the cutoff /x, we still hope that such a dependence would actually be weak. The 
arguments for this case and the corresponding discussion will be presented in the 
last section. 



6 Discussion 

But first let us discuss again the most important question: whether the procedure 
suggested here is really necessary, or perhaps the original procedure of Ref. || is 
still correct and the numerical calculations for the tree-level Hamiltonian (without 
integrating over high momentum modes) provide the results for the correlator ( |1.5| ) 
which do not really depend on ultraviolet cutoff. 

At present, we do not know the answer to this question, and all we can do is 
to provide some arguments both pro and contra regarding the Ambjorn et al. case. 
Let us start with the pro arguments. 

They come mainly from the observation that the physics associated with the 
baryon number non-conservation in electroweak theory at ultrahigh temperatures 
T ^> T c (and the physics of chirality non-conservation in hot QCD) is related to the 
momentum region ~ g 2 T and is essentially non-perturbative. The estimate (|1.7f ) 



18 



for the rate of such non-conservation depends only on the magnetic screening scale 
~ g 2 T and not on the scale gT. It is conceivable then that all relevant momenta in 
the loops for the correlator ( |1.5| ) are of order g 2 T. In that case, the results of the 
classical GR calculation are not sensitive to the ultraviolet region and just do not 
depend on the cutoff \i. 

Essentially the same can be argued in a slightly different way. We have seen that 
the dependence on [i generally comes together with the dependence on the plasmon 
mass scale ~ gT [see Eqs. (|4.6| , |A.19| )1. But the latter is probably not relevant 



because the characteristic momentum scale, where the baryon number nonconserva- 
tion occurs, is of order g 2 T which is much less than m(T) ~ gT. If so, the term ~ 
in m 2 (T,fi) is also not relevant. And that means that the /i - dependent pieces in 
the correlator ( |1.5| ) calculated in a classical procedure are also absent, //this guess 
is true, the results for the correlator ( |1 . 5|) calculated with the tree Lagrangian and 
the results obtained in our procedure with the Braaten-Pisarski-Nair Lagrangian 
should just coincide in spite of that the theories look completely different. The huge 
effects due to effective mass and vertices renormalization should all cancel out in 
the end.0 In this case, all the effective HTL Lagrangian stuff may play the role of 
"(ultra-)violet hands on an enamel wall" |26] - beautiful but irrelevant for physics. 



However, there is also an argument contra. It comes from the analysis of baryon 
nonconservation rate in the temperature region below the phase transition point. In 
this region, the quasiclassical approximation works and an analytical calculation is 
possible ||. When T c -T> gT c , 



r(AB/0) K exp{-^} (6.1) 

where 

m 2 w (T) = m 2 w (0)(l-T 2 /T 2 ) (6.2) 

is the temperature-dependent V^-mass. (At finite temperature, the notion of mass 
is not uniquely defined. The mass (|6.2j) is defined as a curvature of the effective 



8 There are two known examples where the cancellation of the effects related to hard thermal 
loops occurs in a physically observable quantity. One is the axial anomaly [ p7[ . Individual graphs 
for the divergence of anomalous triangle are drastically changed if taking into account the renor- 
malization of Green's function and/or vertices in the heat bath, but these effects completely cancel 
out in the sum of all graphs ( that just follows from the fact that the anomaly d^j^ ~ Tr(GG) is 
an operator identity which does not depend on whether one averages it over vacuum state or over a 
thermal ensemble). A similar cancellation of , at first sight, large effects due to so called anomalous 
damping occurs in the low-frequency electromagnetic polarization operator in quark-gluon plasma 
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potential, not as a pole of the propagator which has a completely different behavior 

Thus, the mass renormalization due to a high frequency thermal loop is essential 
when evaluating T(AB ^ 0) at T < T c . Perhaps, these effects are also essential in 
the region T 3> T c . In this case, the effective HTL Lagrangian is not a luxury but a 
neccessary tool and one should solve the equations of motion with the Hamiltonian 
(|5.8|) to determine the infrared asymptotics of the correlator ( [L.5|) . 

Up to now we neglected the effects due to Higgs and fermion fields. The Higgs 
field is bosonic and can be treated in the same way as the gauge field: one has 
first to calculate hard thermal loops to determine the effective Lagrangian for the 
system including the gauge and the scalar fields and to write it down in a local form 
by introducing when necessary extra dynamic variables. One should then solve the 
equations of motion with this effective Lagrangian for soft modes. 

In fact, Higgs fields have been taken into account in the numerical simulations 
in H. It was found that their inclusion practically does not affect the numerical 
results. That may be thought of as an argument supporting the conjecture that only 
magnetic gauge modes with momenta of order g 2 T are relevant for the problem. 

Fermion fields play a different role because they cannot be treated classically. 
One has to include them however in the loops when deriving the effective Lagrangian 
for the gauge-Higgs system. If the naive GR calculation is ultraviolet-sensitive, also 
the terms in the Lagrangian due to fermion loops are important. If not - they are 
not. 

Our suggestion to people who may wish to determine numerically the behaviour 
of the correlator ( |I.5D at very high temperatures (the calculations in [H] have been 
done in the region T ~ T c ) is the following: First, forget about the Higgs field 
and the fermion fields and do the calculations for a pure Yang-Mills system with 
the ultraviolet cutoff fi ^> gT (but it may be chosen much less than T which is 
a considerable relief for a computer) and look whether the results depend on the 
cutoff. If they do not, the result is probably reliable. However, if they do, one has to 
redo the calculations with the Braaten-Pisarski-Nair effective Lagrangian in the way 
suggested in the paper. If the results obtained with this modified algorithm would 
not depend on the ultraviolet cutoff, one can be sure that they are the answer. 

To find the baryon number nonconservation rate in the hot electoweak theory, 
one has first to derive the analog of Braaten-Pisarski Lagrangian with account of 
Higgs fields and fermions which is yet to be done. 

Finally, let us discuss why one can hope that the modified procedure would lead 
to results which practically do not depend on the cutoff /i. Essentially, it is due to the 
already mentioned fact that /i-dependence in the GR procedure means the relevance 
of high momenta in the loops in which case it is natural to expect that the rate would 
depend not only on the scale ~ g 2 T , but also on the plasmon mass scale ~ gT which 
enters the HTL Hamiltonian ( |5.8| ). And that would mean that the estimate ( |1.7| ) is 
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actually wrong (the heuristic arguments of are physically appealing but they do 
not have the rank of a rigorous proof). In this case, however, one need not bother so 
much with carefullly subtracting the //-dependent counterterms because they affect 
results only slightly (g 2 fiT <C g 2 T 2 ). 

For sure, this argumentation is not rigourous and the possibility remains that 
the estimate (|1.7|) is still correct and the scale ~ gT is irrelevant but the classical GR 
calculation is sensitive to the ultraviolet. In this case one cannot expect that the 
results obtained with the Hamiltonian (|5.8|) are /i-independent, and the only way to 
solve the problem is to accurately isolate /x-dependent counterterms in the effective 
Lagrangian with a gauge-independent (lattice) regularization procedure. We have 
seen that it is a very difficult task, and one should summon the forces to attack this 
problem only if more fortunate options — i) GR procedure is not sensitive to the 
ultraviolet and gives the correct answer and ii) Our modified procedure with the 
Hamiltonian (|5.8|) is not sensitive to the ultraviolet and gives the correct answer - 
are ruled out. 
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Appendix A. Soft gluon vertices with momentum 
cutoff. 

This Appendix is devoted to the calculation of the hard thermal loops for multi- 
gluon vertices with an infrared cutoff for internal loop momentum |p| > fi (gT <C 
fi <C T). Although, as we shall see, the effective Hamiltonian generated by these 
vertices breaks down gauge invariance and cannot be used to solve the problem, 
it is important to understand where the problem is. Also, this calculation may be 
interesting from a methodical viewpoint. We consistently use the real time formalism 
which is simpler and physically more transparent than the imaginary time formalism 
standardly used. 

We start with estimating the contribution from the region |p| > /i in the gluon 
polarization operator. 3 relevant graphs are depicted in Fig. |2]. We work in the 
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simplified Dolan- Jackiw (DJ) version |T7| of the real time technique which is equiv- 
alent to the full-scale real time technique [11]- [14] in this particular problem. The 



D J expressions for the gluon and the ghost propagators have the form (cf . Eq. fll.ijp 




D ahost(P) = ~ S 



(a.i; 



where £ is the gauge fixing parameter. 

The 2-point vertex in the effective Lagrangian is given by the real part of the 
temperature-dependent contribution in the polarization operator which comes from 
the terms where the 5-function insertion in the propagator is taken into account 
only once. In the limit when the external momentum q is considered to be much 
smaller than the internal momentum |p| ~ //, one can neglect g-dependence in the 
vertices and the expression is considerably simplified. Adding 3 graphs together, 
one gets after simple transformations 



n£(g)=20W 



(2vr) 3 (p + q) 2 - p 2 



6[{p + q) 2 ]n B (\p + q \)- 

6(p 2 )n B (\p \)} + O(q) (A.2) 



The imaginary time version of this formula can be found in |Tj|. Note that the 
result is ^-independent. This is the notorious gauge-independence of hard thermal 
loops |L8| . 

Next, we expand in q the denominator (p + q) 2 — p 2 = 2pq + 0(q 2 ) and the 
expression in the square brackets in the integrand in Eq. ( |A.2 ). In the latter, 



there are two terms coming from the expansion of n B (\po + Qo\) and of 5[(p + q) ]. 
When integrating over dpo, we have to take into account two roots of ^-function 
p = ±|p| which give, however, identical contributions. Thus, it suffices to calculate 
the contribution of one of the roots po = |p|. 

The term involving 5'(p 2 ) = ^^^(p 2 ) should be integrated by parts. The 
derivative ^ acts upon ub{po) and also on the structure 



ipyPv - g^p 2 

Po 



(A.3) 



It is not difficult to see that dS flu (p)/dpQ is nonzero only for the spatial components 
H — i,v — j. In the remaining integral over d\p\ the infrared cutoff |p| = e > fi 
should be introduced. The loop integrals for IK 



rah 
L Q0i 



yiab 
11 i0 > 



nff take the form 



yrab 
ii 00 



2g 2 N5 



ab 



1 



dw- 



qv 



g - qv 



00 , 2 dn B 
de e l — - 
de 



(e) 



(A.4) 
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i 

U${q) « 2g 2 N5 ab ^4- I dv I rfee 2 ^ (A.5) 

I 2 J q - qv J M de 




IlfAq) w 2g 2 N5 ab 7 - / dv^vV— ^ / dee 



(2tt) 3 7 [ go - qv 
" ~ J deen B (e)j (A.6) 



with v = p/e. We see that the integrals for n o and Il j factorize in the product 
of two independent integrals: the energy integral and the angular integral. For ELy, 
the expression involves two different energy integrals and the factorization is not 
quite complete. 

Without the infrared cutoff \i the two energy integrals would just coincide - one 
would be obtained from the other by integrating by parts. When \x ^ 0, this is not 
the case: 



- d n B (e) = ^( T 2_3 
de " 3 V 7T 2 



and 



r°° 7r 2 / fi \ 

2 ^ dee n B (e) = — (t 2 - —fjT) (A.8) 

Since these integrals differ by a \i- dependent term, the polarization tensor in Eq. 
(|A.4 ) is not transverse (thereby the gauge invariance is lost). It is rather remark- 
able, however, that this non-transverse piece does not actually depend on the gauge 
parameter £ in the free thermal gluon propagators. 

The appearance of non-transverse piece in the gluon polarization operator with 
momentum cutoff is not a surprise. Just recall that the usual zero-temperature 
photon polarization operator calculated with a momentum ultraviolet cutoff involves 
a quadratically divergent photon mass term. 

What we are going to show next is that the troublesome mismatch ([A. 7 , |A.8| ) 



appears only for the 2-point function. For n > 3 all n-point functions retain the 
same form as without cutoff up to a universal renormalization of the coefficient. 
Consider first the 3-gluon vertex. In Feynman gauge there are two relevant diagrams 
depicted in Fig. |^ (The diagram containing a four gluon vertex does not contribute 
at leading order in q.). Let us calculate the real part of the effective vertex in the 
Dolan-Jackiw technique. Note that here not only the terms involving one thermal 
5-function but also the terms with the product of 3 (^-functions and 3 distribution 
functions contribute. The contribution of these S ® 5 <E> 5 term is huge ~ g 3 T 3 /q 2 har 
where q c har is a characteristic external momentum. It is much greater than the tree 
vertex in the region q c har ~ qT. Moreover, one can see that the integral for this 
contribution is singular in the infrared and is saturated by the internal momenta of 
order q c har- Thus, it is not a hard thermal loop at all! 
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The resolution of this problem is known. The Dolan-Jackiw formalism gives the 
so called T 111 component of the vertex in the matrix real time formalism. What 
we need, however, is the retarded vertex T R (the one which is obtained by analytic 
continuation from the Euclidean region). T R presents the combination of T 111 and 
also other components. The term involving the product of 3 distribution functions 



cancels out in this combination [2J|. One can show that the leading in q terms in 
coincide with those in T 111 which involve only one temperature insertion 0. As 
a result, the integral for the temperature contribution to the 3-point vertex has the 
following structure 



t^a ~ g 3 



d 4 p 

J—^Pl^PluPlX 



S{pl)n B {\p w \) 

^21^31 



ggjWjPgoj) S(p 2 3 )n B (\p 30 \) 

012032 



13O23 



(A.9) 



Pk ~Pi an d Phi fe = 1, 2, 3 are the internal momenta in the loop. 5ki have 
Tq c har ~ gT 2 and are small. We have also neglected the dependence 



where Ski 

the order ~ 1 q char 

on the external momenta in the tensor structure multiplying the square bracket in 



The individual terms in the r.h.s. of Eq. (|A.9 ) have the order ~ g 3 T 3 /q 2 har and 



are large. However, these large contributions cancel in the total sum due to the 
algebraic identity 



53(1,2,3) 



1 



+ 



1 



21031 



012032 



+ 



1 



0~13°~23 







(A.10) 



As a result, the vertex has the normal HTL order ~ g 3 T 2 /q c h ar . To get this "non- 
leading" term, one has to expand the numerators of the fractions in Eq. ( |A.9|) in the 
external momenta with respect to, say, the momentum p\. As earlier, there will 
be the terms coming from the expansion of the distribution functions /ib(|P2o|) and 
^b(|P3o|) which are proportional to external energies and also, potentially, the terms 
coming from the expansion of 5-functions. The terms ~ 8'(pi) are dangerous as, 
when integrating by parts, the derivative may act upon the structure ~ P^PuPx/po 
and provide the contribution involving the integral ( |A.8| ) rather than ( |A.7| ). That 
could give rise to some extra structures in the vertices and, as a result, in the effective 
equations of motion( TA.18| ). 

Fortunately, it does not happen. Really, the terms ~ S'(p 2 ) in ( |A.9|) are multi- 
plied by 



1 

$32 



1 

0~23 







(All) 



Now we understand what almost in the discussion preceding Eq.(1.9) meant 
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Thus, only the terms involving the expansion of n,B{\p2o\), n B{\P3o\) and thereby 
proportional to external energies and involving only the integral ( |A.7|) appear. 

Let us prove it in general (for a similar proof in imaginary time framework see 
PHI). The proof is inductive. Suppose we already know that, for the n-point vertex, 
the large individual contributions ~ g n T 3 / (g c har) n_1 coming from the graphs where 
5-function term is inserted in a particular internal line cancel in the total sum, 
and terms coming from the expansion of the 5-function also cancel out. Let us prove 
it for the (n + l)-point vertex. 

The relevant graphs are shown in Fig. [?]. The integral has the form 



r n+i f d4 P 



(27T) 



5{pj)n B (\pio\) 

021 ■ ■ ■ O n +l,l 
^(Pn+l) ra B(|Pn+l,o|) " 
°l,n+l • • • °~n,n+l 



(A.12) 



Let us prove that 



S n+1 (l,...,n + l) = - + - + ••• 

021 •• ■ On+1,1 °12 ■ ■ ■ n +l,2 

+ - — = (A.13) 

°l,n+l • • • °n,n+l 

To this end, present the second term in the r.h.s. of Eq.( |A.13|) as 

11 1 

+ 1 H x (A. 14) 



0\2 ■ ■ ■ °n+l,2 ^12^32 ■ ■ ■ O n +l,l 0~l,n+\0~32 ■ ■ ■ n +l,2 

and do the same for the third, . . ., and n-th term in the sum. Then it is not difficult 
to see that 

. . . , n + 1) = — !— [S n (l, . . . , n) - S n (2, . . . , n + 1)] (A.15) 

°n+l,l 

which is zero by the inductive assumption. 

To get a nonzero result, one has to expand o~(Pk)n>B(\pko\) in the integrand in 
([A. 12|) in the external momenta. The term proportional to 5'(pf) involves the factor 

' 1 -S n (2,...,n + 1) = (A.16) 



0~32 ■ ■ ■ <Wi,2 $2,n+l ■ ■ ■ S n ,n+1 

Thus, an additional non-transverse structure involving the difference of the integrals 
(|A.8|) and (|A.7|) appears only in the gluon polarization operator. Unfortunately, is 
it already bad enough to spoil the game. 

10 As earlier, only the terms with a single temperature insertion contribute in the retarded (n+ 1) 

- point vertex. 
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Let us try to write down the effective equations of motion corresponding to the 
n-point functions calculated above. By means of the time contour formalism of 
Sect. 3 one can see that indeed the retarded Greens functions enter the equations 
of motion ensuring causality. Thus we obtain 

D.F^ix) = f(x) (A.17) 



D^{x) = f(x) + ^N^TA\x) (A. 18) 

The current j M is given by eqs. ( p.6| ), ( |5.7| ) with the plasmon mass m p i(T) replaced 
by m p i(T,ii), where 

mJ(T,/i) = \g 2 N (t 2 - ^//t) (A.19) 



The extra term oc A 1 in eq. ( |A.18| ) appears due to the mismatch ( |A.7| )), ( |A.8| ) 



discussed above. Due to this term the equations of motion ( |A.1S ) are not consistent 



with Gauss' constraint ( |A.17| ). This means that choosing initial conditions satisfying 



the constraint (|5.11|) , one cannot assure that this constraint will be fulfilled also at 
later times, and the classical problem cannot be consistently posed. 



Appendix B. Thermal loops on a Hamiltonian lat- 
tice. 

We present here some illustrative calculations which display the way the cutoff- 
dependence arises in the classical GR procedure with the lattice ultraviolet regular- 
ization. Unfortunately, the results of this study are negative — for the gauge theory 
in interest we were not able to extract explicitly the cutoff-dependent terms which is 
necessary to construct an explicit lattice-based cutoff-independent algorithm. But 
it is important to understand where the problem is. It is conceivable that after 
considerable future efforts such an algorithm would be eventually constructed. 
Consider first the \(f) A theory. The lattice Lagrangian is 

4- = ££#-? E (0n +ei -0n) 2 -^ 3 E^n (B.l) 
Z n A n;i=l,2,3 4- n 

n It is more convenient to calculate graphs in the Lagrangian formulation, but one should re- 
member that our lattice is 3-dimensional and the time is real and continuous 
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where ei = (1, 0, 0) etc. To calculate graphs, we have to go over into the momentum 
representation (see for details). Actually, the only thing we need to know here 
is the modified dispersive law 

P 2 = Po " 4 S 1 - ™s(a Pi )] = (B.2) 
a i=i 

The 3-momenta p, range within the limits 

— n/a < pi < n/a (B.3) 

Let us calculate the graph in Fig. [I] on the lattice in the real time formalism. The 
real time scalar lattice propagator has the form (cf. Eq.( |1.9|) ) 



In the limit aT ^> 1, we can, as earlier, take the classical limit T/\p \ of the Bose dis- 
tribution function after which the thermal contribution to the scalar mass acquires 
the form 



' 4 J-n/aJ-n/aJ-n/a (2tt) 3 [3 - ^Ll COS ( a Pi)] 



With the help of the identities 



1 r°° 1 

dae~ aA , — I dxe acosx = J„(a), (B.6) 



A Jo 2vr 

the integral can be reduced to a one-dimensional form as spelled out in Eq. 

Let us go over now to the gauge theories. For illustrative purposes, we restrict 
ourselves with the case of scalar QED. The part of the lattice Lagrangian involving 
scalar fields is 

£ = a 3 ]T|X? o n | 2 + a Y, (CfeA+e*, n0n + 
n n;i=l,2,3 

<Pn +Sl U: +ei , - 0; +ei 0n +ei - <P*M (B.7) 

where V = d — ieA 0n and U n+ei>n are complex fields living on the links of the 
lattice . (They have the meaning of parallel transporters: t/ n+ei , n ~~ > exp{ieaAj(x)} 
in the continuum limit.) 

The temperature contribution to the photon polarization operator is determined 
by the graphs in Fig. ^. In the continuum limit, it is given by the same integral as in 
Eq. (|A.2|) . Let us impose the lattice ultraviolet cutoff and assume that aT ^> 1 ^> aq 
where q is the external momentum. 
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Consider first the component Hoo of the polarization operator. It is a little bit 
simpler than others because the vertices are the same as in the continuum theory. 
For n o the only modification brought about by the lattice are the limited range of 
integration ( |B.3|) over spatial momenta and the modified dispersive law ( P-2| ). As a 
result, the integral for n 00 reads 

r^/ a dp 2 



if 



00 



2e 2 T / djh 



■*l a dpi 

-7r/o 2lX 



*7 a dp 3 2p\ - p 2 



Tv/a 2n (p + q) 2 — p 2 
'S[(p + q) 2 ] 5{p 2 ~ 



(B.8) 



bo + go I |po| 

(where as earlier we neglected small q compared to p ~ 1/a whenever possible). 
Expanding l/\po + go I an d 5[(p + q) 2 ] in g, we get for the term oc 1/a 

2e 2 T r r r d 3 s Y,i gisinsj 1 



n 



00 



where Sj = api and 



U! 



f (2vr) 3 q - ^- a J2i Qismsi u 3 
2 -2^(1 -cos Si ) 



(B.9) 



And now we see that n o is not the function of only q and q 2 = q\ + q\ + q\ as is the 
case in the continuum theory, but depends on each component of q separately. It is 
more or less obvious from the form of the integral. But if a reader cherishes a hope 
that it might not be true, he is welcome to expand the expression ( |B.9|) over q/|go| 
up to the fourth order and to be explicitly convinced that on top of the rotationally 
invariant structure ~ (q 2 + q 2 + g|) 2 also the structure ~ q\q\ + q\q\ + qfq 2 pops 
out with a non-zero coefficient. 

For completeness, we present here also the results for the other components of 

n, 



n 



2e 2 T 



0i 



-go 



dh 



sin Si 



TT (2tt) 3 g - rr£i gisins^ 



1 



IF 



2e 2 T 



-go 



dh 



sin Sj sin s 



1 



, (2tt)3 g 



i 



Efc gfcsinsfc uj 



(B.10) 



(B.ll) 



As the lattice regularization does not break gauge invariance, Yl^ u of Eqs. (p.9| 



B.ll|) is transverse: g^Il^^ = 0. 



But, as the rotational invariance is broken, the cutoff-dependent counterterm in 
the 2-point vertex of the effective HTF Fagrangian is not presented as the sum of just 
two standard (transverse and longitudinal) tensor structures as was the case in the 
effective theory without cutoff. The same refers to counterterms for multiple-gluon 
vertices which we have not calculated. 
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cf. the expression (|A.4|) with the momentum cutoff. 
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Figure Captions 



Figure 1 The contribution to the propagator which yields the effective temperature 
dependent mass in A0 4 theory. 

Figure 2 The polarization contribution to the vector boson propagator at finite tem- 
perature. The dashed lines represent ghost propagators. 

Figure 3 The sum of self energy insertions which yield the propagator in \<p A theory. 

Figure 4 The contour of the time integral in the path integral representation for the 
correlation function. 

Figure 5 Two-loop ladder graph contributing to the effective action for very soft mo- 
menta q ~ g 2 T. 

Figure 6 Leading contributions to the effective three gluon vertex in Feynman gauge. 
The dashed lines represent ghost propagators. 

Figure 7 Hard thermal loop for a (n + 1) - gluon vertex. The vertical dashes stand for 
the temperature (5-function) insertions. There are also similar graphs with 
internal ghost loop. 

Figure 8 One-loop contributions to the photon polarization operator in scalar QED. 
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